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Abstract 



O 



We study in this letter the double beta decay of ^^^^Xe with emission of two neutrinos which has been recently measured 
^by the EXO-200 collaboration. We use the same shell model framework, valence space, and effective interaction that 
■we have already employed in our calculation of the nuclear matrix element (NME) of its neutrinoless double beta decay. 
'Using the quenching factor of the Gamow- Teller operator which is needed to reproduce the very recent high resolution 
.^■^^Xe (^He, t)^'^^Cs data, we obtain a nuclear matrix element M^''=0.025 MeV~^ compared with the experimental value 
;M2''=0.019(2) MeV-i. 
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The double beta decay is a rare process (second or- 
■der in the weak interaction) which takes place when the 
[single beta decay of the parent even-even nucleus to the 
■neighbor odd-odd nucleus is forbidden by energy conser- 
[vation or highly suppressed by the angular momentum se- 
■lection rules. In addition it is one of the major sources of 
^background for the even rarer neutrinoless decay which, if 
■detected, will settle the nature (Majorana or Dirac) and 
[the mass scale of the neutrinos. Until the EXO-200 mea- 
■sure [ij -'^'^^Xe was the only (experimentally relevant) po- 
[tential neutrinoless emitter whose two neutrino decay was 
■unknown. In addition, the lower bound to its half life pub- 
lished by the Dama collaboration ^] demanded a nearly 
■complete cancellation of the nuclear matrix element. Af- 
[ter the EXO measure, we know that the matrix element is 
■small (indeed, the smallest among the measures ones) but 
not pathologically so (see in Table [T] the present status of 
■the 2v decays from the recent compilation of ref. Q). The 
[EXO-200 measure has been confirmed by KamLAND-Zen 
■ 0] only a few weeks ago. 

The 2v decay half-life contains a phase space factor and 
the square of a nuclear matrix element 



The nuclear structure information is contained in the 
nuclear matrix element to which only the Gamow- Teller 
at^ part contributes in the long wavelength approxima- 
tion. 
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(0;^|(Tt+|m)(TO|CTi+|0+) 
- (M, + Af/)/2 
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Therefore, to calculate the nuclear matrix element we 
need to describe properly the ground state of the par- 



Table 1: Experimental 2v fifi decay matrix elements 



Decay 



M^^-^) (MeV- 



48Ca^ 


48Ti 


0.047±0.003 


4.4 X lO^^ 


^^Ge ^ 


™Se 


0.140±0.005 


1.5 X 10^1 


82Se-> 


82Kr 


0.098±0.004 


9.2 X 


9621- _^ 


96Mo 


0.096±0.004 


2.3 X 1019 


looMo ^ iooRu 


0.246±0.007 


7.1 X 10^8 


"6Cd - 


> "6sn 


0.136±0.005 


2.8 X 101*3 


i28Te _ 


► i28xe 


0.049±0.006 


1.9 X 1024 


i30Te - 


► i30xe 


0.034±0.003 


6.8 X 10^0 


i36Xe - 


J. i36Ba 


0.019±0.002 


2.1 X 10^1 


isoNd - 




0.063±0.003 


8.2 X IQi*^ 



ent and grand daughter nuclei as well as all the 1+ ex- 
cited states of the intermediate odd-odd nucleus. In other 
words, the GT~ strength function of the parent, the GT+ 
strength function of the grand daughter and the relative 
phases of the contributions from each intermediate state. 

Our description of the wave functions of the states 
involved in the process is based in the Interacting Shell 
Model approach. The valence space includes the orbits 
0g7/2, ld5/2, ld3/2, 2si/2, and Ohn/a, covering the sector 
of the nuclear chart between N,Z=50 and N,Z=82. We 
use the effective interaction gcn50:82 [5|] which is based 
in a renormalized G-matrix obtained from the Bonn-C Q 
potential using the methods of ref. Q • The final interac- 
tion is obtained through a (mainly monopole) fit to about 
300 energy levels from ~90 nuclei in the region with a rms 
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deviation of 100 keV. More details can be found in ref. [8[. 

It is well known that the effective Gamow- Teller oper- 
ator a for complete harmonic oscillator valence spaces 
can be approximated hy q ■ a t^. q is called the quench- 
ing factor and behaves as a sort of effective GT charge (see 
ref. 9] for a recent update of this topic). The value of q has 
been fitted throughout the nuclear chart and the resulting 
values are 0.82, 0.77, and 0.74 for the p, sd and pf shells. 
Asymptotically it tends to 0.7. With the quoted value of 
q for the pf-shell the half-life of the 2i/ double beta decay 
of ''^Ca [10| could be predicted in perfect agreement with 
the later measured value fll'|. The problem arises when 
we try to describe heavier emitters in which the minimal 
complete valence spaces (in the harmonic oscillator sense) 
are still out of reach computationally. A possible solu- 
tion is to carry out a fit to all the experimentally available 
GT decays. We did this exercise in our valence space and 
the resulting value was q—0.57. A more accurate way of 
estimating the quenching factor is by comparing the theo- 
retical predictions for the Gamow Teller strength functions 
relevant for the process with the experimental results ob- 
tained in charge exchange reactions. These data were not 
available for the "'^^^Xe case until the advent of the re- 
sults of the ^'^^Xe ("^He, t)-'^'^^Cs reaction in the appropriate 
kinematics, which have been published very recently [l2j . 
These results impact in our calculations in two ways; first 
because they give us the excitation energy of the first l"*" 
state in ^^^Cs, 0.59 MeV, unknown till now, which appears 
in the energy denominator of equation [51 and secondly be- 
cause it makes it possible to extract directly the quenching 
factor adequate for this process. 



tions in the same valence space and with the same effective 
interaction which we had used in our calculation of the 
matrix elements of ^^'^Sn, ^^STe, "OTe and ^^^Xe in ref. 
13[. First, we compare in Fig. [T] the theoretical running 
sum of the B(GT~) strength of ^'^^Xe with the experimen- 
tal data from [12] • We have normalized the total theo- 
retical strength in the experimental energy window to the 
measured one. This implies a quenching q=0.45. Notice 
the very good agreement between the theoretical and ex- 
perimental strength functions. If we had shifted the the- 
oretical position of the first excited state of ^■^^Xe to its 
experimental value, the quenching factor would have been 
slightly larger. 
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Figure 1: (color online) The running sum of the Gamow- Teller 
strength of -'•^^Xe (energies in MeV). The theoretical strength is nor- 
malized to the experimental one. 



In what follows, we use for the A=136 isobars the wave 
functions that result of the large scale shell model calcula- 



Figure 2: The running sum of the 2^ matrix element of the double 
beta decay of ^^^Xe (energies in MeV). 



Then we compute the 2^ matrix element with the 
quenching factor extracted above. The result is given in 
Figure [H in the form of a running sum. The final matrix 
element M^'^=0.025 MeV~^ agrees nicely with the experi- 
mental value. However, one should bear in mind that the 
absolute normalization of the Gamow- Teller strength ex- 
tracted from the charge exchange reactions may be affected 
by systematic errors, which could lead to modifications 
of the extracted quenching factor. Minor variants of the 
gcn50:82 interaction which locally improve the quadrupole 
properties of "^Ba lead to q=0.48 and M2''=0.021 MeV'^ 
Comparing Figure [T] and Figure [5] it is evident that the 
nuclear matrix element M^" does not saturate in the ex- 
perimentally studied energy window. In fact, according to 
the calculation, about 40% of the total 2^ matrix element 
comes from states above it. 

For completeness, we present in Table [5] a compilation 
of the 21/ matrix element elements for which there are large 
scale shell model calculations. For the results of QRPA- 
like calculations see refs. [13, In the calcula- 
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Table 2: The ISM predictions for the matrix element of several 2v 
double beta decays (in MeV~^). See text for the definitions of the 
valence spaces and interactions. 







IVl l^tAp J 


q 


M^'^^'tb'l 
iVl I Lll 1 


TNT 






0.047±0.003 


0.74 


0.047 


kb3 




48rpi 


0.047±0.003 


0.74 


0.048 


kb3g 




48Ti 


0.047±0.003 


0.74 


0.065 


gxpfl 


™Ge^ 


76 Se 


0.140±0.005 


0.60 


0.116 


gcn28:50 


™Ge^ 


76 Se 


0.140±0.005 


0.60 


0.120 


jun45 




82Kr 


0.098±0.004 


0.60 


0.126 


gcn28:50 




82Kr 


0.098±0.004 


0.60 


0.124 


jun45 




i28xe 


0.049±0.006 


0.57 


0.059 


gcn50:82 


i30Te ^ 


i30Xe 


0.034±0.003 


0.57 


0.043 


gcn50:82 


"6Xe ^ 


i36Ba 


0.019±0.002 


0.45 


0.025 


gcn50:82 



tion of the decay of ^®Ca, the valence space is the full pf 
shell and the interactions are defined in refs. [3] (kb3), 
[111 (kb3gV and ^ (gxpfl). The details can be found in 
refs. [H, 21 1 . The quenching factor conies from a fit to all 
the experimentally available Gamow- Teller decays in the 
region using the kb3 interaction [g^]- No equivalent fits 
are available for kb3g or gxpfl, therefore we use the same 
quenching factor for all of them. We are convinced that 
the differences would be negligible in this case. 

For the decays of 76Ge and ®^Se we take a core of ^6^1 
and the valence space spanned by the orbits lp3/2, Ofa/ii 
lPi/2j and Ogfl/ ^. T he interactions are defined in refs. [23| 
(gcn28:50) and [24] (jun45). There is a published calcula- 
tion with a preliminary version of jun45 in ref. (25| . The 
value of the quenching factor q=0.60 was obtained in this 
reference from a fit to the Gamow- Teller decays in the 
region. As there are no fits available with the other in- 
teractions we adopt this value for them all. The results 
for the 21^ NME compare quite well with the experimental 
values and are very much interaction independent, as in 
the ''^Ca decay. 

In addition, we have plotted in Figure [3] the running 
sum of the matrix element of the 2v double beta decay 
of ''^Ca, using the kb3 interaction, and in Figure |3] the 
same for the decays of 76 Ge and *^Se using the gcn28:50 
interaction. Notice that the patterns in the three case are 
quite different and are also at variance with the results for 
the ^^6xe decay shown in Figure [21 reflecting the different 
structures of the intermediate nuclei. The ^^Ca decay is 
peculiar in one aspect, to be the only case in which there 
are important canceling contributions to the 2^ NME. In 
fact, if the contributions from all the states would have 
had the same sign, the matrix element would have more 
than doubled, thus needing a quenching factor q=0.53 in- 
stead of q=0.74, in line with our findings in the other va- 
lence spaces. In the remaining cases the sign coherence is 



nearly complete. However, it is difficult to make a more 
precise surmise about the extra quenching factor needed 
when incomplete valence spaces (in the harmonic oscillator 
sense) are used, based solely in this example. Concerning 
the energy at which the 21/ NME saturates, the differences 
among the four cases studied are relatively minor, and one 
can safely conclude that states beyond 10 MeV excitation 
energy above the first 1+ state in the intermediate odd-odd 
nucleus do not contribute to M^"^. 
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Figure 3: The running sum of the 2u matrix element of the double 
beta decay of **Ca (energies in MeV). 
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Figure 4: (color online) The running sum of the 2u matrix element 
of the double beta decays of ^^Ge (red) and *^Se (blue) (energies in 
MeV). 



Finally, for the decays of ^^^Te and ^^^Te we use the 
same valence space and interaction than for the decay of 
^'^6Xe. The quenching factor, as we have already men- 
tioned, is obtained by a fit to the single Gamow- Teller 
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decays experimentally known in this region (and within 
our computational limits). Had we used the newer value 
from the analysis of the -'^'^^Xe charge exchange data, the 
values of the matrix elements would have been: 

M2'^(i28Te)=0.037 MeY'^ 
M2'^(i30Te)=0.027 MeV"! 

In conclusion, we have shown that large scale shell 
model calculations which describe in detail the spectro- 
scopic properties of large regions of the nuclear chart can 
also make accurate predictions (or postdictions) of the 
nuclear matrix elements of weak processes including the 
rarest ones; the 2i/ and Ov double beta decays. In the 2i/ 
case we have found that we can explain the experimental 
data provided the Gamow- Teller operator is renormalized 
(quenched) so as to reproduce the Gamow Teller single 
beta decays or the charge exchange results in the relevant 
regions. We have highlighted the case of ^^^Xe, whose 2^ 
double beta decay half-life has been recently measured by 
the EXO-200 collaboration and confirmed by KamLAND- 
Zen. 
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